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We consider non-linear generalizations of fractal interpolating functions applied to functions of 
one and two variables. The use of such interpolating functions in resizing images is illustrated. 
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I. INTRODUCTION 

An Iterated Function System (IFS) may be used to construct fractal interpolation functions for some data 
The simplest example of interpolating a function x(t), given data points (U, Xi), i = 0, 1, . . . , N, starts with an IFS 
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The coefficients a n , c n , e„, and f n determined from the conditions, for n = 1, 2, . . . , N, 
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which leads to 
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| With these, the transformation of Eq. |QJ can be written as 
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Wn[t ) = i> = t n + ^4 tn-l 
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W„(x) = x' = — f — r x n + — — - X n -1 (4) 

^ ' \tn ~ t-n— lj (,%— 1 — t n ) 

in which form it is apparent W n {x) = x' is determined by a linear (in t) interpolating function between the points 
I_h ' (in-ii ^n-i) and (t n , ar„). Graphs of fractal interpolating functions can then be made by applying the random iteration 
5^ , algorithm: 

• initialize (t, x) to a point in the interval of interest 

• for a set number of iterations 

- randomly select a transformation W n (t, x) 

- plot {t',x') = W n (t,x) 

- set (t,x) = (t',x') 

• end for 
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In this paper we consider two non-linear generalizations of such fractal interpolating functions. The first concerns 
how to extend the linear interpolation of Eq. to higher-degree interpolations. The second generalization arises 
when one considers the construction of fractal interpolating functions for functions of two (or more) variables - in 
this case, even a linear interpolation of the form of Eq. fllf. when applied to each variable, will result in a non-linear 
interpolating function. This has an obvious application in representing two-dimensional images, where a function of 
two variables (the pixel coordinates) could represent a black-and-white image (using a Boolean function), a gray-scale 
image (using a scalar function), or a colour image (using a vector-valued function of the three rob Jred, green, blue] 
values) . This problem has been examined extensively in the context of image compression 0, 0, S j m the last 
section we consider a related problem of using these iterated function systems to rescale images, or portions thereof. 



II. FUNCTIONS OF TWO VARIABLES 



We first consider a function z(x, y) of two variables, and examine the problem of constructing a fractal interpolating 
function from the data Xj, yj, Z{ j, where i = 0, 1, ... , M, j = 0, 1, . . . , N, and Zi j = z{xi, yj). To this end, consider 
the transformation 

Wmn(x) 

^mn(|/) — CmnU fmn 

W mn (z) = A mn x -f~ B mn y C rnn xy -h D 7nn (5) 
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We then impose, for m = 1, 2, . . . , M and n = 1,2, ... ,N the conditions 
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The coefficients turn out to be 
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With these, the transformation of Eq. © can be written as 
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in which form it is apparent T / F mn (z) = is determined by a function implementing a linear interpolation over the 
grid (x m -i 5 y n -i), (a? m -i,y n ), (x m ,y n _i), and (x m ,y n ). 
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III. QUADRATIC INTERPOLATING FUNCTIONS 

The interpolating functions considered up to now have used a linear interpolating formula between adjacent points 
to construct the IFS. In this section we indicate how this can be generalized to quadratic interpolations. 

A. Functions of one variable 

For a function x(t) of one variable, using data points (ij, x{), i = 0, 1, . . . , N, consider the transformations 

W n (t) = a n t + e n 

W n (x) = c n t + d n t 2 + f n (10) 
and impose the conditions, for n = 2, 3, . . . , N, 

W n ( t0 ) = ( tn ~ 2 ), W n ( tM ) = ( tn - 1 ), W n ( tN ) = ( tn ). (11) 
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The point tu is determined as 
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with corresponding point x m . The coefficients of the IFS are determined as 
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with which the transformation can be written as 
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In this form we see that a quadratic (in t') interpolating function is used between the points (t n ,x n ), (t n -i,x n -i), 
and (t n -2,x n - 2 ). 

B. Functions of two variables 

We next consider a function z(x,y) of two variables, and construct a fractal interpolating function which employs 
a quadratic interpolation between points. To this end, consider the transformation 
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We then impose, for m = 2, 3, . . . , M and n = 2, 3, . . . , N, the conditions 
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along with the corresponding z points. The coefficients of the IFS can then be determined, by which the transformation 
of Eq. Hlfty can be written as 
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Although tedious to work out, the generalization of the preceding considerations to higher-order interpolating func- 
tions is straightforward in principle. 



IV. IMAGE SCALING 



As an application of the preceding, in this section we consider the task of rescaling a colour image. This is a natural 
problem for an interpolating function of two variables (x, y) interpreted as pixel coordinates - the function z(x, y) in 
this case will be a vector-valued function having three components representing the rgb value of the pixel specifying 
the amount of red, green, and blue present. 

The procedure used to scale an image of size M pixels wide by N pixels high is as follows. We first read in the 
rgb values of each pixel of the image, and use that as the data to construct a fractal interpolating function z(i, j), 
where i = 1, 2, . . . , M and j = 1, 2, . . . , N. To then resize the image, so that the resulting image is of size s x M x s y N, 
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we construct a new fractal interpolating function z'(i,j) = z(s x i, s y j). Applying the random iteration algorithm to 
choosing independently a transformation index at each stage, will then result in the rescaled image. The 
generalization of this procedure to rescale a portion of an image is straightforward. 

As examples of the results of this procedure, consider the figures in the Appendix. We start with the image 
appearing in Fig. ^ and zoom in on the area of the face. The result appears in Fig. [21 together with a comparison 
done using a simple linear interpolation scheme. Zooming further into the area of the eye results in Fig.|3 again with 
a comparison of the result of a simple linear interpolation. Generally, the number of iterations needed in the random 
iteration algorithm to produce acceptable images is of the order of s x M x s y N, where the original image is of size 
M x N. Also, while slower, the quadratic fractal interpolating function typically produces, for the same number of 
iterations, a "smoother" looking image than the corresponding linear interpolating function. However, as with all 
interpolation schemes, there comes a point where such higher order interpolating formulas actually start to produce 
worse results due to an artificially high sensitivity to fluctuations in the data. 

Some informal tests of this procedure seems to indicate that better results are obtained for images of people, natural 
scenery, etc., as opposed to those containing lettering, simple geometric shapes, and similar constructs. This might 
be expected, given the general fractal nature of such objects in nature. However, as with all interpolating functions, 
it is important to remember that no structural information beyond that of the original image is being provided (for 
example, one could not zoom in on the face of Fig. ^ to such a degree as to see individual skin pores). 

The preceding demonstrates that these non-linear fractal interpolating functions of two variables can be used in 
principle to represent images. It would be interesting to extend these considerations to the case of partitioned iterated 
function systems, upon which much work has been done with respect to compressing images \4$. Work along these 
directions is in progress. 
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Scaling of figures 




FIG. 3: 



Enlargement of the area around the eye of Fig. Q via a) a fractal interpolation function b) a linear interpolation 
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